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1 Introduction 

The aim of the LHCb experiment at the Large Hadron Colhder at CERN is to perform 
precise tests of the Standard Model (SM) in the flavour sector, in order to disentangle 
possible New Physics (NP) effects. In particular the measurement of the branching 
fractions (BF) of the rare leptonic decays of the B, D and K mesons can give hints 
of the presence of NP particles at the tree and loop levels. The LHCb detector is well 
suited to study decays with muons in the final states: very efficient trigger allows to 
collect events containing one or two muons with very low transverse momenta; very 
good momentum resolution 5p/p = (0.4-0.6)% reflects into excellent invariant mass 
resolution {a{M) ~ 25 MeV for B two-body decays); the offline muon identification 
permits to have a good muon efficiency, e ~ 90%, for a muon misidentification rate 
less than 1% for 1 < p < 100 GeV/c. 

In this paper the analyses of four purely leptonic rare decays with ~ Ifb^^ of pp 
collisions collected by LHCb in 2011 at ^/s = 7 TeV are presented. 

The rare decay Ks — )■ f^'^f^^ are a very useful source of information on the short- 
distance structure of AS* = 1 Flavour Changing Neutral Current (FCNC) transi- 
tions. This decay is suppressed in the SM and the prediction on its BF is B{Ks — )■ 
fi~^fi~)sA4 = (5.0 ± 1.5) X 10^^^ P [2]. The current best limit on this decay, obtained 
in 1973, is equal to B{Ks /i+At") < 3.2 x 10"^ at 90% of C.L. [3]. The contributions 
of NP to the BF, e.g. from light scalars, are allowed up to one order of magnitude 
above the SM expectation. 

A blind analysis is performed on 1 fb~ of data collected during 2011 that contains 
~ 10^'^ Ks inside the LHCb acceptance. The main sources of background are due 
to combinatorial muons from semileptonic decays and to Ks — > vr+vr" where both 
pions are misidentified as muons. The Kl — )■ fj^~^fi~ component is negligible for 



this analysis. The LHCb mass resolution is exploited to discriminate the Ks — )■ 
TT+TT^ with both pions misidentified as muons. Moreover, to increase the signal and 
background separation a multivariate classifier, a boosted decision tree (BDT), based 
on geometrical and kinematic informations is used. The number of expected signal 
events, for a given branching ratio hypothesis, is evaluated by normalizing to the 
Ks — )■ 7r+7r~ events. This normalization reduces the common systematic uncertainties 
between the two channels. The modified frequentist method, CL^, is used for the 
upper limit determination The CLs curves for B{Ks f^^fJ'~) are shown in 
Figjl| The observed upper limit is: B{Ks < 9.0 x 10"^ at 90% C.L., with 
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Figure 1: Expected CLs (dashed blue line) under the hypothesis to observe 
background-only. The green (dark) band covers 68% (la) of the CLs curves obtained 
in the background-only pseudo-experiments, while the yellow (light) band covers 95% 
{2a). The solid line corresponds to the observed CLg. 

an improvement of a factor ~ 30 with respect the previous best limit. 

3 B^^ ^ 

The search for B^^^^ — )■ is one of the best way for the LHCb experiment to 

constrain the parameters of NP models with an extended Higgs sector. These decays 
are highly suppressed in the SM because they are both FCNC and helicity suppressed. 
The SM model predictions are [5J: i3(5° ^ yU+/i~)5M = (3.2 ± 0.2) x lO'^ and 
B{B^ — 7- fJ'^fi~)sM = (1.0 ±0.1) X 10^^". If we consider the Minimal Supersymmetric 
SM (MSSM) in the high tan/3 approximation, the enhancement of B{B^ — >■ n'^fi~) 
due to scalar and pseudoscalar sector is found to be proportional to tan^/3 [B]. 

The results presented here are based on ~ lfb~^ of data recorded by LHCb in 
2011 [7j. The events are classified using a BDT, based on kinematic and geometrical 
variables, and the dimuon invariant mass. In order to avoid any bias, the mass region 
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Figure 2: Expected CLg (dashed black line) under the hypothesis to observe 
background-only (left) for the — )• jj,^ ji~ and CLs under background-plus-signal 
events according to the SM rate (right) for 5° — > /U^/U^, with yellow area covering 
the region of ±lcr of compatible observations; the observed CLg is given by the blue 
dotted line; the expected (observed) upper limits at 90% and 95% C.L. are also shown 
as dashed and solid grey (red) lines. 



'rr^^J,^l = [m{B^) — 60 MeV/c^, m(i?°) + 60 MeV/c^] is blinded until the analysis is final- 
ized. The number of expected signal events, for a given branching ratio hypothesis, is 
evaluated by normalizing to the decays — )■ Kit, B^ — )■ J/ipK^ and 5^ — )■ J/ip(j). 
The compatibility of the observed events with background-only and background plus 
signal hypothesis is then computed. The modified frequentist method CLs is used 
for the upper limit determination The CLs curves for B{B^ — )■ fi~^fj^~) and 
B{Bg — )■ /i"^/x~) are shown in Fig|2] 

The upper limit obtained are: -6(5° /x+yU") < 4.5x 10~^, and B{B° yU"*"/^") < 
1.03 X 10^^ both evaluated at 95% C.L.. In order to compare the upper limit on 
B{B'^ — )■ with the theoretical prediction, this value has to be multiplied by 

0.911 ± 0.014, which takes into account the effective lifetime of the Bs meson [8j. 

These two upper limits have been also recently combined |TTj with the results 
obtained by CMS with 4.9 fb"^ [9] and ATLAS with 2.9 fb"^ [lOj. The combined 
upper limits are: S(5° /U+/i") < 4.2 x 10"^ and B{B^ /W+At") < 8.1 x 10"^° 
both at 95% of C.L.. Despite the fact that the i3(i?° — )■ /i^/U~) limit is very close 
to the SM prediction, NP contributions are still possible. For example, destructive 
interference between NP and the SM could decrease significantly the value of the 
B{B^s^^^^^^-). 

The — ^~ branching fraction is very suppressed in the SM. Its main contribu- 
tion is due to two-photons intermediate state which sets the SM prediction bounds: 
10"^^ < i3 < 6 X 10"^^ at 90% C.L. [12]. Many NP models predict enhancements for 
this BF, e.g. RPV-SUSY, which predicts B ~ 10~^ due to a tree level transition [13]. 
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The current best experimental limit is i3 < 1.4 x 10~^ at 90% C.L., from Belle [H]. 

An analysis has been performed on a data sample of 0.9 fb~^ and selecting the 
decay chain D*^ — )■ D^{-^ fi'^fi~)7c^. The selection of the events is based on geometric 
and kinematic properties of and its daughter particles. After this selection the 
background is mainly composed of combinatorial muons from semileptonic decays of 
b and c hadrons and the peaking background of D*^ — )■ Z)°(— )■ h~^h~)7r^ events, with 
double hadron misidentification. The combinatorial background is reduced making 
use of a multivariate discriminant based on geometrical and kinematical information. 
The signal events are normalized to the D*^ — )■ )■ 7r"'"7r~)7r^ channel. Since 
the same selection is applied to signal and normalization channel, this procedure 
permits to reduce common systematic uncertainties in both channels. The event 
yield is determined from a 2D fit on the dimuon invariant mass and the difference 
between the D*~^ mass and mass. A preliminary result is given: i3(Z)° — t- < 

1.3 X 10^ at 95% C.L. [B]. 

5 B^^-^ — > fi~ jj,^ 

The 5°^) — )■ jj,^ jjL^ jj,^ is a FCNC process and takes its largest contribution from the 
resonant decay — > J /ipcf) in which both mesons decay into two muons. In the SM 
a non-resonant process can also occur through the exchange of a virtual photon with 
a BF ~ 10~^° [IE]. However, NP effects might enhance the B{B^g-^ — t- fi^ fi" fi^) by 
the exchange of new particles at the tree level. 

A cut-based selection has been designed on the 1 fb~^ data sample collected during 
the year 2011. The selection algorithm is tuned using the resonant decay mode 
and optimized keeping blind the signal region. The combinatorial background was 
estimated from the mass sidebands. The selection criteria for signal and control 
channels are based on particle identification, separation between the B vertex and 
the primary vertex, the quality of the B decay vertex. A veto on and B^ masses 
is applied in order to exclude B~^ — )• J/ipK^ and B — )■ K*(f). All the non-resonant 
peaking background yields in the signal region are found to be negligible. 

The signal BF is measured by normalizing to B^ — )■ J/iIj{^ fi^ fi~)K*^{^ K^ix^) 
decays selected with the same criteria. Systematic uncertainties are evaluated by 
comparing Monte Carlo simulation with data. 

After unblinding, one event is observed in the B^ signal window and no events are 
observed in the Bg window. The number of observed events is consistent with the 
expected background yields. The CLg method has been used to evaluate the upper 
limits on BF [17|: B{Bl^^ fi+ fx- fi+ fx-) < 1.3 x 10"* and B{B^ -> fx+ fx" fx+ fx-) < 

5.4 X 10-9 both at 95% C.L. 
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6 Conclusion 



We presented some of the purely leptonic rare decays analyses based on 1 fb~ col- 
lected by LHCb during the 2011 run. A brand new result on the upper limit on 
B{Ks — )■ /i^/U") is presented. This result improves the previous best limit by a factor 
~ 30. The results on the upper limits on B{B^ — )• /i^/i^) and -23(5° l^^l^ ) Put 
severe constraints on NP and further improvements are foreseen. The combination 
of the upper limits obtained by ATLAS, CMS and LHCb has been also reported. 
This combination improves on the limits obtained by the individual experiments and 
represents the best existing limits on these decays. The study on the upper limit 
on B{D^ —7- /i^yU~) shows an improvement on the previous best limit by an order 
of magnitude. A search for the decays 5° — ?■ fi^fi~fi^fi^ and — )■ ^^fi^^^fi~ is 
also discussed and the upper limits on their branching fractions are reported. Many 
improvements are foreseen in the near future. With the 2012 run LHCb expects to 
more than double the data set and will put even more stringent constraints on the 
phase space of many NP models. 
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